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SUMMARY 


A Gaussian multilayer diffusion model is used to calculate nine independent 
solid rocket motor (SRM) exhaust cloud dispersion cases in order to characterize 
potential environmental impact. These cases are based on seven standard meteorologi- 
cal regimes for the Cape Canaveral, Florida, area representing characteristic over- 
land advection cases. The dispersion cases include the following: A modified "spring 

fair weather" Space Shuttle case, with about 60 metric tons HCl exhausted below 2 , 5-km 
altitude (a conservatively large HCl source strength) ; all seven meteorological cases 
for the chemically similar but smaller Titan III vehicle, with resultant HCl source 
strengths up to 15 metric tons for stabilized clouds bounded at 2,0 km; and an 
abnormal (pad abort) Titan III case. Downwind vertical HCl concentration profiles 
are analyzed, and vertical HCl column density O is found to be well characterized 
by the decay expression O = aX”^ when the downwind distance from launch site X is 
defined as the sum of the downwind distance from the cloud stabilization point X^g 
and an effective virtual source distance upwind from the cloud stabilization 
point Xq. (a and 3 are empirical constants,) The dispersive O decays are 
illustrated, and sets of a, X^, and resultant HCl source strength are tabulated 

for use with either X or elapsed time t. Equivalent decays of vertically averaged 
HCl concentration C(HCl) are also presented. 

The calculated decays of O and C(HCl) differ greatly among the seven meteo- 
rological regimes, A range of more than two orders of magnitude in O and C (HCl) 
is spanned at X > 100 km and t ^ 2.0 hr. At shorter distances, the total span 
in a and C(HCl) still exceeds an order of magnitude for all X > 10 km and 
t > 0.2 hr. These results suggest that various meteorological conditions at Cape 
Canaveral lead to widely different exhaust cloud dispersion rates. Also, values of 
a ^ 6000 ppmv-m and C(HCl) ^ 5 ppmv are calculated 'for X - 50 km and t - 5.0 hr 
for the two least dispersive Titan III cases. This could result in acidic rainwater 
of pH ^ 1.5, which can cause significant damage to some plants. 

The set of calculated C(HCl) decays is compared with published analyses of 
in-cloud peak HCl concentration data from eight Titan III launches. Four major fea- 
tures are noted. First, straight-line fits of the log C(HCl) versus log t data 
are in agreement with the calculated power-law characteristics. Second, the data 
exhibit a similar spread (range) of two orders of magnitude for t > 1.0 hr, depending 
on launch meteorology. Third, an inclusive envelope that bounds the calculated decays 
of C(HCl) also bounds nearly all the in-cloud HCl concentration data for 
t ^ 0.2 hr. Finally, the measured decays of in-cloud HCl tend to be less dispersive 
than the calculated ones. This apparent difference is consistent with the use of 
model variances based on relatively small-scale (compared with SRM cloud size) turbu- 
lence measurements. Despite this last difference, we conclude that the calculated 
data agree well with the experimental data. In effect, this provides a partial vali- 
dation of the model as an assessment tool, since individual model-prediction versus 
experimental-measurement comparisons are not made. 

Several factors affect the validity of the calculated O and C(HCl) values, 
especially at large X and t. First, the model assumes stable stratification con- 
ditions in the lower troposphere and, hence, lack of vertical convective motion. 

This oversimplification invalidates use of the model under many conditions. Second, 
neglect of convective loss of HCl from the SRM exhaust cloud’s upper boundary, HCl 



sorption at ground level, and variable advection for the respective MDM layers tends 
to result in unrealistically large values of O and C(HCl) at large X. Neglect 
of these factors, however, is a convenient means of allowing the model to predict 
somewhat high at large X, so as to err on the safe side. Finally, the above- 
mentioned use of small-scale variances tends to result in underestimation of O and 
C(HCl) at large X. Evaluation of this effect is beyond the scope of this paper . 

Since the basic formulation of the modified Gaussian dispersion models is 
statistical-analytical in nature and does not treat the essential physics of cloud- 
scale and mesoscale dynamics, a need is apparent for resuming development of a much 
more comprehensive model- Such a predictive tool is clearly needed for assessment of 
atmospheric effects, especially in the vicinity of the Cape Canaveral land-sea inter- 
face where sharp contrasts in surface thermal and moisture characteristics exist and 
significant convective activity occurs. 


INTRODUCTION 

The National Aeronautics and' Space Administration (NASA) has been examining the 
possible environmental impacts of its Space Shuttle Program for about 8 years. Formal 
Environmental Impact Statements were published in July 1972 (ref. 1) and April 1978 
(ref. 2) . One of the potential problem areas cited, atmospheric pollution, stems 
mainly from the planned use of a solid rocket motor (SRM) booster design. The tropo- 
spheric portion of this environmental problem centers on the possible effects of 
relatively large, localized, low-level releases of SRM exhaust products. For the 
present discussion, these include more than 60 metric tons of hydrogen chloride (HCl) 
and about 90 metric tons of aluminum oxide (alumina) particles emitted below 2 . 5-km 
altitude per launch.^ 

Environmental impact studies of SRM exhaust clouds in the troposphere have been 
shaped and focused by three principal concerns. First, gaseous components of the 
afterburned exhaust (refs. 3 to 5) combined with chloridized alumina particles 
(refs. 6 to 13), entrained ground debris (refs, 14 to 16), and large amounts of 
sprayed water used for acoustic baffling and launch pad cooling (refs. 5 and 14), may 
adversely impact ground receivers, which include both plant and animal life (refs. 1, 
2, and 17 to 21) . Second, the possibility exists that precipitation scavenging of HCl 
might lead to localized deposition of unacceptably acidic (mostly hydrochloric acid) 
rain on nearby land areas or protected waters before atmospheric dispersion has 
effectively reduced the potential hazard to acceptable limits (refs. 1, 2, 6, and 
17 to 24) . Third, the potential exists for inadvertent weather modification effects, 
which conceivably may occur up to a few days after launch (refs. 25 and 26) and may 
be caused by aliomina ice nuclei (refs- 13 and 25 to 29). 

The formation, altitude stabilization, and transport with atmospheric dispersion 
of Space Shuttle and Titan III SRM exhaust clouds have been the subjects of several 


^Note that the stabilized ground cloud is defined in ref. 2 as containing 
35 metric tons HCl (20 sec burn up to about 1.1-km altitude), and more recent mission 
designs (1976) indicate about 44 metric tons HCl exhausted for a 20 sec burn up to 
1.2-km altitude. Since the ground cloud definition in ref. 2 was used for evaluating 
surface level concentrations, additional HCl exhausted at levels higher than 1.2 km 
would be relatively unimportant. However, for the present Shuttle case (prediction 
of in-cloud concentrations and eventual acid-rain characteristics) , a more conserva- 
tive resultant source strength (61 metric tons HCl) is being used. 
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modeling studies (refs. 14 and 30 to 35). Experimental field studies of Titan III 
clouds (refs. 1 , 8, 10 to 13, 27 to 29, and 36 to 48) have attempted to provide cloud 
geometry and temporal species concentration data for validation and refinement of 
NASA's Gaussian multilayer diffusion models (refs. 30 to 35) and also for more 
sophisticated atmospheric dynamics models currently under development (refs. 22, 23, 
and 49) . The field studies, in conjunction with idealized laboratory experiments 
(refs. 6, 7, and 9), have also provided certain physical-chemical data to allow 
initial characterizations of cloud composition and important rate processes for calcu- 
lation of acid-rain formation and weather-modification effects. Although the early 
Titan III field studies concentrated on obtaining transient-species concentration 
data at ground level from arrays of stationary and semimobile sites (refs. 36 to 41), 
more recent efforts have focused mainly on obtaining in-cloud species concentration 
profiles and particle/aqueous-acid aerosol size distribution measurements from suc- 
cessive figure-eight patterns flown by instrumented aircraft (refs. 39 to 48). 

Although most of the published modeling results have focused on defining 
transient-species concentrations and dosages at ground level, relatively few have 
dealt with corresponding species concentrations aloft, even though vertical concen- 
tration profiles are inherent in the calculations and in-cloud concentrations have 
been measured. Since studies of potential acidic rain and inadvertent weather- 
modification effects require information on HCl concentrations (and alumina acid- 
aerosol characteristics) aloft, an initial objective of this study is to characterize 
representative sets of downwind vertical HCl concentration profiles that apply to 
large exhaust clouds over a range of meteorological conditions. 

For certain applications, it is both practical and desirable to define integrated 
vertical burdens, or vertical coliamn densities, of components. For example, optical 
remote sensing techniques (two-ended or Sun tracking) , which lead directly to measure- 
ment of column density, have a practical advantage in mapping the source strength 
characteristics of a polluted cloud since only two-dimensional data in the horizontal 
plane are required. Thus, calculations of both vertical concentration profiles and 
column densities may be useful in the consideration of certain remote sensing 
applications . 

At present, however, information about column density in SRM exhaust clouds is 
needed to calculate acid-rain deposition characteristics. It has been shown, through 
development of an idealized washout model (refs. 6 and 24), that scavenging of HCl (g) 
by raindrops is effectively an irreversible absorption process that is linearly 
dependent on vertical HCl (g) column density C, regardless of the actual variation in 
vertical concentration profile below the rainfall source. Thus, useful analytical 
expressions for predicting HCl(g) washout, rain pH, and ground deposition of HCl can 
be applied to independently dispersing SRM exhaust clouds if the temporal and spatial 
variation of a can be specified in terms of time and ground coordinates. 

Accordingly, the primary focus of the present paper is to evaluate vertical HCl 
column density under a variety of meteorological conditions for input to the refined 
acid-rain model of reference 24. Some pertinent details of the acid-rain model are 
first reviewed, however, and then some specifics of the present approach are given. 

The analytic rain scavenging model treats the idealized case of an independently 
generated vertical rainfall that overrides and scavenges HCl (g) , by washout processes, 
from an independently dispersing SRM exhaust cloud. The washout model applies best 
at low-to-moderate relative humidities for SRM exhaust clouds, where HCl (g) tends to 
predominate over the aqueous-acid aerosol component after a few minutes of cloud 
dilution, and at stable stratification conditions in the lower troposphere. The 
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washout model was first used (ref. 6) to predict the resultant characteristics of 
acid-rain deposition for a "spring fair weather" (SFW) Space Shuttle SRM exhaust 
cloud dispersion case, derived independently from application of Model 4, one of the 
Gaussian multilayer diffusion models (MDM) that were developed at the NASA Marshall 
Space Flight Center (MSFC) . The SFW dispersion case was based on one of seven 
"standard meteorologies" documented for the Cape Canaveral, Florida, area. These 
were selected to represent a range of characteristic cases of overland advection with 
turbulent diffusion in the planetary boundary layer in order to establish an initial 
basis for assessing the downwind concentration history of HCl at the Earth's surface 
(refs. 50 and 51). 

In the present paper, the Model 4 MDM is applied to nine independent exhaust 
cloud dispersion cases based on the seven standard meteorological regimes. These 
cases include the following: (a) An improved version of the SFW Shuttle case; 

(b) all seven standard meteorological cases for the smaller but chemically similar 
Titan III SRM-propelled vehicle, which exhausts about 40 percent of the material that 
the Shuttle does in the 0- to 4.0-km altitude range; and (c) an abnormal (pad-abort) 
Titan III case. Downwind vertical HCl concentration profiles calculated for these 
cases provide one basis for comparisons with existing in-cloud data on Titan III 
launch effluents. They also provide the requisite data for determinations of a as 
described below. 

A simple approach is used herein for quantitative characterization of O from 
extensive arrays of multilayer HCl concentration histories calculated with the MDM. 

It is then shown that each dispersive decay of O can be generalized as a function of 
downwind distance (or time) from launch site through a one-term power-law expression. 
Sets of two empirically derived constants for each of the nine dispersion cases are 
shown to characterize the power-law decays of O within ±10 percent for the entire 
range of MDM output distances downwind of SRM exhaust cloud stabilization (1 to 
100 km). The O expressions, in terms of the previously developed HCl(g) washout 
model, are defined as follows (refs. 6 and 24): 

r ^top _o 

^ " I " Pq(HC1)z - ax ^ (1) 

Jo 

where p^ ^(HCl) is the prewashout HCl(g) concentration in parts per million by 
volume at^altitude z, Pq(HCI) represents a vertically averaged concentration 
over z, X is the downwind distance from launch site, and a and 3 represent 
empirical constants. These O expressions are sufficiently simple, when used as 
inputs for the analytic acid-rain model, that they allow relatively straightforward 
parametric analyses and characterizations of downwind acidic rainfall for assumed 
occurrences of rainfall events. 


SYMBOLS 

C(HCl) calculated average HCl (g + aq) concentration in diluted SRM exhaust cloud, 
expressed in dimensionless volume/volume units, ppmv 

/\ 

C(HCl) peak value of measured in-cloud HCl concentration, expressed in dimension- 
less volume/volume units, ppmv 

m^ initial mass of HCl in SRM exhaust cloud, g HCl 
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p(HCl) 


average concentration of HCl (g) in parcel of diluted SRM exhaust, expressed 
in dimensionless volume/volume units, ppmv 

p^(HCl) initial value of p(HCl) before HCl (g) washout begins at t - 0, ppmv 

Qj effective heat release for solid rocket propellant combustion, cal/g 

propellant (1 cal = 4.184 J) 

t elapsed time after launch, hr 

mean unidirectional wind speed which describes horizontal motion of SRM 
exhaust cloud, m/sec 

X unidirectional downwind distance from launch site, + ^cs' 

X^ virtual source distance upwind from SRM cloud stabilization point at 

X^^g = 0, km 

X^g unidirectional downwind distance from SRM cloud stabilization point, km 

SRM cloud-centroid height, m 
height of SRM cloud top above ground, m 

a, B constants which define power-law decay of HCl (g) column density O in 

dispersing SRM exhaust cloud; a is equivalent to ppmv-m at X = 1 km, 
and B is dimensionless exponent of X in equation (1) 

0^ average wind direction, deg 

A washout coefficient for HCl scavenging by falling raindrops, sec“^ 

a vertical HCl (g) column density, ppmv-m 

Abbreviations : 

CFP cold front passage 

FFW fall fair weather 

FW, Pre-CF fair weather, pre-cold front 

LLSB low-level sea breeze 

MDM multilayer diffusion model 

Post-CFP post-cold front passage 

SB sea breeze 

SFW spring fair weather 

SRM solid rocket motor 
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APPLICATION OF THE MULTILAYER DIFFUSION MODEL 


The rocket effluent dispersion calculations used in this paper employ the 
NASA MSEC exhaust cloud rise preprocessor model and the Gaussian multilayer diffusion 
models (MDM) , which are described and documented in references 30 and 31 for opera- 
tional prediction of toxic fuel hazards. Model 4 (multilayer, diamond-shaped stabi- 
lized cloud) , Version II (updated enthalpy content) of the MDM code is used through- 
out this study for dispersion calculations and is designated herein as MDM-4(II). 
Simultaneous inclusion of the precipitation-scavenging subroutine is indicated by the 
designation MDM-5(II) (ref- 30). 


Model Description 

An overview which broadly characterizes the NASA MSFC rocket exhaust diffusion 
modeling technique is given in section IV. A of reference 33. Specific applications 
of the cloud rise preprocessor and the variously modified Gaussian multilayer dif- 
fusion models to rocket exhaust dispersion for tropospheric air quality predictions 
at ground level are described in references 14 and 30 to 35. Some of these refer- 
ences focus on the inherent mathematical and physical limitations of this basically 
statistical analytic approach (ref. 14) , parametric studies of the principal variables 
upon which the MDM*s are based (refs. 14 and 32), and comparisons of calculated rocket 
effluent concentrations at ground level with those obtained from other advective dif- 
fusion models (ref. 14) . 


Standard Meteorological Regimes for Cape Canaveral 

The meteorological regimes selected for study consist of the set of standard 
meteorologies (ref. 50) that were originally provided for calculating atmospheric 
dispersion of rocket exhaust effluents in the Cape Canaveral, Florida, area (ref. 51) . 
They represent the major types of meteorological conditions which lead to overland 
transport and are likely to be encountered in the Cape Canaveral area. 

The meteorological profiles of temperature, wind speed, and wind direction from 
references 30 and 51 for all seven meteorological regimes are shown in figures 1 (a) 
to (g) . They are based on averages of selected data obtained from Kennedy Space 
Center (KSC) rawinsonde releases and from the NASA 150-Meter Ground Winds Tower 
Facility at KSC. The profile data, including pressures, are tabulated in 
reference 50. 

Inspection of the meteorological profiles indicates large differences in the 
thermal structure and wind characteristics aloft. Furthermore, note that significant 
deviations occur, at various levels, from the preprocessor-calculated average values 
(one each) of wind speed U^ and wind direction that are used in MDM-4(II) to 

describe average SRM-cloud advective motion after altitude stabilization. Thus, 
horizontal wind shear effects are not directly accounted for in the MDM-4(II) calcu- 
lations since individual layers in the multilayer structure are not allowed to advect 
independently, even though average wind speed values for each layer were specified in 


9 

The "KSC Sea Breeze, Normal Launch" case in ref. 50 is the present low-level 
sea breeze case. The present sea breeze case is not tabulated in ref. 50, but is 
illustrated in ref. 30. 
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the program output. Instead, the effects of horizontal and vertical wind shear, 
relative to average SRM exhaust cloud advective speed and direction, are simulated in 
each layer by sets of horizontal (cross direction) and vertical turbulence parameters 
which applied to each layer. These parameters are standard deviations of the wind 
azimuth and elevation angle fluctuations. In principle, these standard deviations 
can be obtained from a comprehensive set of local meteorological data, such as those 
tabulated in part for the Cape Canaveral area on specific launch occasions (refs. 52 
to 54) . In the absence of direct measurements, the turbulence parameters can be 
deduced from the simple profile measurements such as those shown in figures 1(a) 
to (g) . The latter procedure was used for the seven standard meteorologies used in 
this study. The methodology used to deduce these parameters is described in refer- 
ence 14 and references 30 to 35- The specific values of MDM-4(II) input parameters 
used in this study are summarized in table A1 in the appendix . 


Characteristics of Stabilized SRM Source Clouds 

The SRM exhaust cloud rise preprocessor numerically calculates the vertical and 
horizontal distribution of HCl which results from buoyant cloud rise with turbulent 
air entrainment under a specified potential temperature gradient (refs. 30 to 35). 

The vertical source-strength distribution at cloud altitude stabilization is assigned 
a multilayer structure, based partly on locations of changes in the meteorological 
profiles. Similarly, the stabilized cloud is assigned a specific geometric shape in 
the vertical cross section, which characterizes the layer widths, and hence the HCl 
concentrations, in the stabilized SRM cloud. In the case of Model 4, a conical body 
of revolution, having a diamond- shaped cross section symmetrical about the vertical 
centroid axis, is used to bound the horizontal layers. This is shown in figure 2 and 
is discussed below. 

An example of the geometric definition of a Titan III stabilized SRM exhaust 
cloud is shown in figure 2(a) for the sea breeze (SB) meteorological regime at Cape 
Canaveral- Resultant layer source strengths are also shown for completeness. The 
calculated height of the cloud centroid is 993 m and is based on a currently 

accepted value for the effective propellant heat release Qj of 2790 cal/g propellant. 
The fine-layered definition of the surface mixing layer has an upper bound of 

800 m based on the onset of temperature inversion at 800 m as shown in figure 1(d). 

The "stem" of the so-called column begins at 1680 m and has a radius of 200 m. 

A second stabilized SRM exhaust cloud is shown in figure 2(b) for the Post-CFP 
meteorological regime applied to the case of a normal launch of the Titan III vehicle. 
In this case, the calculated height of the cloud centroid is 1341 m, compared with 
1506 m for the "pad abort" launch case (Post-CFP (pad abort) ) , and 993 m for the 
previously discussed SB case. Thus, as expected, both the vertical temperature pro- 
file and the input source strength influence the determination of Z^. Note also that 
the height of the surface mixing layer (H^^^ = 1400 m, from fig. 1(g)) is greater than 
in the SB case but is not defined strictly at the minimum temperature point (1700 m) 
in figure Kg). Instead, is defined at the turning point in potential 

temperature . 

A third Titan III SRM exhaust cloud is shown in figure 2(c) for the CFP case. 

The corresponding vertical temperature gradient in figure 1(f) is not as steep as in 
the Post-CFP case, resulting in a somewhat lower cloud centroid height (Z^^ = 1230 m) . 
However, the temperature profile does not exhibit an inversion up to 2000 m, and thus 
the surface layer is assigned a height of H^^ = 2000 m; the multilayer structure is 
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arbitrarily terminated at this point- Geometric definitions for the remaining meteo- 
rological cases (not shown) are based on application of the same principles as 
described in the above three cases. Plots of the layer structures and initial verti- 
cal distribution of HCl for these remaining cases (Oj = 2790 cal/g propellant) are 
similar, based on tabulated values of HCl source strength in each SRM exhaust cloud 
layer Q(HCl) in the appendix and an entrainment coefficient Y of 0.64 for 
stabilized Titan III exhaust clouds. 


Source-Strength Considerations 

Although the vertical-line — source-strength input function (propellant mass burn 
rate versus altitude) is accurately known for any prescribed Titan III or Space 
Shuttle mission, the dynamics of cloud rise determine the resultant vertical and hori- 
zontal distributions of SRM exhaust mass (and concentrations) in the altitude- 
stabilized cloud. Thus, the calculated layer source strengths, tabulated in the 
appendix for the present MDM- subdivided stabilized SRM exhaust clouds, differ from 
the incremental exhaust input source strengths and are highly dependent on the meteo- 
rology considered. Consequently, the resultant cumulative HCl source strength is 
always smaller than the cumulative input of HCl mass exhausted up to various alti- 
tudes- Also, both the resultant layer and cumulative HCl source strengths vary sig- 
nificantly for different potential temperature profiles. Illustrative examples of HCl 
input and resultant source strengths are discussed in the remainder of this section. 

A summary of the resultant cumulative HCl source strengths is given in table 1; other 
quantities in this table are discussed in the next two sections. 

A graphical summary of HCl input histories and resultant source strengths for 
the chemically similar SRM boosters for both the Space Shuttle (Nov. 1973 design) and 
the Titan III launch vehicles is shown in figure 3(a); figure 3(b) shows the inputs 
for two more recent Shuttle missions (1976 design) . The cumulative HCl mass inputs 
to the MDM for both the Shuttle and the Titan III cases are based on propellant con- 
sumption rates and launch mission trajectories. They differ by a factor of 2.44 up 
to 2.0-km altitude in figure 3(a). The resultant "preprocessed" fall fair weather 
(FFW) cumulative HCl source strength in the MDM stabilized cloud for Titan III is 
also shown in figure 3(a) for altitudes up to 2.0 km, where the calculation was 
terminated- Resultant layer source strengths are tabulated for all the Titan III 
cases in the appendix. 

The original Space Shuttle spring fair weather (SFW) cumulative source strength 
for the stabilized cloud, used in previous acid-rain calculations (ref. 6) , is shown 
as a single point in figure 3(a). In view of the Mission 2 November 1973 design 
(ref. 54) and also the more recent 1976 design (refs. 55 and 56) shown in figure 3(b), 
it should be reduced to about 60 metric tons of HCl- This figure approximates cumu- 
lative HCl exhausted up to 2 . 5 km and is still consistent with retention of the 
original, more conservative 4.0-km upper cloud boundary for the Space Shuttle (as 
opposed to 2.0-km altitude bounds used for the Titan III cases). 


Vertical HCl Profiles and Downwind Cloud Dispersion 

The maximum centerline HCl concentration for each defined MDM cloud layer is 
selected from the MDM-4(II) output in order to provide conservative upper limit HCl 
concentrations for the respective layers, and thus for the vertically integrated col- 
umn densities. Sets of calculated HCl concentrations thus define vertical profiles 
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for nine sequentially increasing distances downwind of the SRM cloud stabilization 
point ranging from 1 to 100 km. 

Vertical HCl concentration profiles at six values of are shown in fig- 

ure 4(a) for the FFW case. The plotted HCl concentrations, which apply mostly at (or 
very near) the respective centroids of each defined layer, are arbitrarily connected 
by straight-line segments in this semilog plot- For ^cs - ^ most of the HCl is 
contained within layers 6 to 10, and the vertical concentration profile in layers 1 
to 5 quickly approaches uniformity in this active and well-defined surface mixing 
layer. A zero— absorption (perfect reflection) surface boundary condition is used 
throughout this study. The concentration just above layer 10 is shown to be 
0.001 ppmv. In actuality, this concentration is zero, since the layers are termi- 
nated at 2000 m to define the upper cloud boundary. This boundary is also designated 
as a perfect reflector. Thus, no mechanism for convective loss at the boundary of 
the top layer or sorption loss at the surface has been provided. 

Graphical representations of the downwind vertical HCl concentration profiles 
for the remaining seven meteorological cases are shown in figures 4(b) to (h) , using 
the same parametric range for ^cs figure 4(a). Inspection of these various 

cases indicates large differences in initial HCl concentration profile at X^g = 1, 
in subsequent downwind decay of concentration within respective layers, in downwind 
development of a uniform surface mixing layer, and in development of upper cloud 
layer stratifications which sometimes are relatively nondispersive with respect to 
decay of HCl concentration. The development of characteristic surface mixing layers 
can be seen to vary in height from 200 to 2000 m and in the rate of approach to 
vertically uniform HCl concentration, although uniformity always occurs within 
Xcs ^ 20. 

For the Titan III launches under the seven standard meteorological regimes at 
Cape Canaveral and for the pad abort case, table 1 summarizes the SRM exhaust cloud 
dimensional, advection, altitude-stabilization, and source-strength characteristics 
for the present MDM-4(II) application. Note that the mean wind speed for the cloud, 
which correctly relates distance to time in the model, differs from the wind speeds 
of the individual layers (tabulated in the appendix) . The cloud stabilization times, 
shown from the model output, are used to estimate downwind drift distances from the 
launch pad to the point where cloud altitude stabilization theoretically occurs. 

These calculated cloud stabilization drift distances are compared later with empiri- 
cally determined virtual source distances. 

Finally, as indicated earlier, the resultant HCl source strengths for the 
stabilized SRM exhaust clouds are also summarized in table 1. These apply from the 
Earth's surface to the MDM upper cloud boundary, and are conserved throughout the SRM 
cloud dispersion history (in the absence of rain) , since no loss terms are applied 
at the SRM cloud boundaries . 

APPLICATION OF CALCULATED HCl PROFILES 
Determination of Vertical HCl Column Density 

Since the fundamental HCl source parameter in the idealized acid-rain model of 
references 6 and 24 is vertical column density O, the MDM-4(II) cloud dispersion 
results can be incorporated into the acid-rain model but still remain uncoupled, pro- 
vided that an independently integrated expression for O can be determined. This 
effectively allows replacement of a large set of vertically dependent (multilayer) 
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concentration terms in a three-dimensional acid-rain model with a single expression 
independent of the vertical coordinate variable. 

Initially in this study, single-term power-law expressions for O were postu- 
lated as the simplest conceivable cases, even though significant curvature was found 
(ref. 6) for the previous SFW Space Shuttle case. These were explicitly dependent on 
downwind horizontal distance of the cloud centroid from the cloud stabilization point. 
The problem, therefore, is to determine whether or not the proposed power-law 
expressions for O adequately characterize the calculated results for the various 
meteorological cases in figure 4. 

Integrated vertical HCl column densities were obtained as follows. Tabulated 
vertical profiles of peak centerline HCl concentration, identical to those shown in 
figure 4 but including three additional downwind distances were integrated 

numerically at each by a product summation technique. The resultant summations 

were numerically equivalent to the area under a linear concentration versus altitude 
curve obtained by connecting adjacent-layer centerline concentrations with straight- 
line segments. The C values were initially plotted on a log-log scale as a function 
of downwind distance from cloud stabilization (X^g from 1 to 100 km) . Significant 
curvatures were found at X^g 5 20 km for each of the meteorological cases, indi- 
cating departures from simple power-law behavior similar to that found in the previ- 
ously analyzed SFW Space Shuttle case (ref. 6) . 

In an effort to find a simpler and more useful characterization of O that 
would still retain the necessary accuracy, it became evident that addition of a 
relatively small empirically determined distance X^ to X^g led to a good straight- 
line fit of log a versus log (X^^g + Xq) over the entire range of X^g for all 
the meteorological cases studied. Since all the empirically determined Xq values 
were positive and reasonably small, falling entirely in the range 2.5 to 5.0 km with 
a single exception (FFW) requiring 10.0 km, X^ was identified as an effective 
virtual source distance for cloud stabilization. These results for X^ indicate 
reasonably good (but not perfect) correspondence with independent estimates of SRM 
cloud drift between vehicle launch and the MDM-4(II) specified cloud stabilization 
times. (See table 1.) Accordingly, the corrected distance of the cloud from launch 
site is defined by X = X^^g + Xq for the present study. 


Results of Vertical HCl Column Density Determinations 

The empirical power-law fits of O versus X are shown in figures 5 (a) to (h) 
for each of the Titan III cases studied. In addition, the previous (ref. 6) and 
presently modified SFW Space Shuttle cases are shown in figures 6 and 7. The corre- 
sponding decay expressions O = a{X^^ + Xq)”^ are given in each case, along with the 
resultant HCl source strength m^ and the mean transport wind speed U^, which cor- 
rectly relates downwind distance to elapsed time. 

The straight-line fits of the calculated data points in figures 5(a) to (h) and 7 
are excellent throughout the data range for all the meteorological cases studied. 

Thus, the demonstration of adequate power-law fits for downwind decay of 0 leads to 
substantial analytic simplification, with respect to the calculation and parameteriza- 
tion of acid-rain characteristics, through application of the previously defined 
analytic expressions. 
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In order to examine and compare some of the resultant effects of meteorology 
on O, the above-mentioned empirically determined parameters are summarized in 
table 2 along with other relevant exhaust cloud properties derived from the MDM-4(II) 
calculations. The corresponding time-dependent expressions for O, calculated 
through use of to convert distance to time, are also given. 

Although the deduced values of Ot can be referred to as an HCl column density 
” source strength" at unit distance, the reader should recognize that a is inherently 
sensitive to Xq. In turn, is influenced in the MDM- 4 ( I I ) cloud rise calcula- 

tions by the vertical profiles of temperature and wind speed, as well as the place- 
ment of vertical layer boundaries- Thus, the physical significance of a is not 
precisely defined. 

A noteworthy internal consistency exists for the two Post-CFP meteorological 
cases, which apply to a normal launch (line source) and a pad abort (point source) 
situation. The respective Xq and 3 values are identical, and the ratio of a 
and a values (0.71) is identical to the ratio of the m^ values- In the absence 
of other similar paired calculations for Titan III launches, it is not known whether 
a simple linear scaling law would apply to other meteorological cases. 

A roughly similar comparison can be made between the fully modified Space 
Shuttle (see fig. 7) and Titan III SFW cases, but it is not as exact since additional 
multilayer structure bounded at 4000-m altitude is used for calculating the Shuttle 
case, whereas the Titan III case is terminated at 2000 m. While the respective 
values of 3 differ (1.64 and 1.98), the Xq values are the same (4.0 km), and the 
ratio of values of a (3.85) is fairly close to the ratio of values of mQ (4.21). 

In terms of corrected distance from launch site, the point at which the ratios of O' 
and mQ are equal to 4.21 is X = 1.3 km, which is still close to the launch site. 
Thus, the above comparisons between Titan III and Space Shuttle launches for SFW 
meteorology suggest that approximately linear scaling of O with source strength 
applied. 


Collective Results for Dispersive Decay of O and Potential Rain pH 

The power-law decays of O, deduced in figures 5(a) to (h) for the eight 
Titan III cases, are shown collectively as functions of corrected distance from 
launch site X in figure 8 and also as corresponding functions of elapsed time after 
launch t = X/(3.6U^) in figure 9. The SFW Shuttle case deduced in figure 7 is 
shown as a dashed line in each figure. Rain pH values shown in figures 8 and 9 were 
calculated from reference 24 and are discussed later. 

The most noticeable feature of figures 8 and 9 is that the dispersive decays 
of O differ greatly among the seven standard meteorological regimes- A range of 
more than two orders of magnitude in O is spanned at X - 100 km downwind and/or 
t ^ 2.0 hr. At shorter distances and times the total span in O is somewhat less, 
but still exceeds an order of magnitude for X > 10 km and/or t > 0.2 hr. These 
results tend to confirm our earlier expectations of large variations in atmospheric 
dispersion under widely different meteorological conditions, and they help to empha- 
size the need for developing a comprehensive atmospheric dynamics model (e.g., see 
refs. 22 and 49) in order to deal effectively with this large source of variability. 
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A second important feature of figures 8 and 9 is that relatively large O values 
and low potential rain pH*s are shown to be possible out to relatively large downwind 
distances and correspondingly long elapsed times. For example, a pH - 1.5 is con- 
sidered environmentally significant for single exposures (refs. 19 and 20) . If we 
accept the validity of the two meteorological cases which tend to define upper bounds 
for O (FFW and FW, Pre-CF) , initial rain pH*s between 1.0 and 1.5 (depending on 
rainfall rate) could occur at downwind distances up to 50 km and elapsed times of 
nearly 5 hours, based on O' ^ 6000 ppmv-m (ref. 24) . Although these estimates are 
derived from the two least-dispersive meteorological cases, they are not considered 
excessively conservative. For example, somewhat more severe stagnation conditions 
are possible, and no estimated uncertainty bounds have been imposed. Although the 
predictions in figures 8 and 9 also indicate O > 2000 ppmv-m and potential rain pH*s 
of 2.0 or less at downwind distances up to 200 km and elapsed times exceeding 10 hr, 
we must recognize that the validity of these longer range results becomes less certain 
with increasing distance and time. 

Some competing effects, which both reduce and increase O, are now identified . 
First, since convective loss of HCl from the SRM exhaust cloud's upper boundary, HCl 
sorption at ground level, and variable advection (wind speed and direction) for the 
respective MDM cloud layers are not accounted for in the present model calculations, 
inclusion of these processes would tend to reduce O' systematically at progressively 
large distances. Second, the variances used in the present MDM-4 (II) calculations 
were originally based on field data which applied to relatively small-scale atmo- 
spheric turbulence (ref. 51) . It is well known (ref. 57) , however, that as the rele- 
vant turbulence scale increases for atmospheric dispersion of very large plumes, the 
appropriate variances and overall dispersion rates tend to become weaker power-law 
functions of the characteristic scale size (e.g., source cloud diameter) and downwind 
distance. Thus, the exponent 3/ which characterizes the decay of a with X, could 
systematically decrease to about 0.50 at large values of X under some meteorological 
conditions. This would occur in a fashion similar to that observed in large-scale 
dispersion studies (ref. 57) and that actually observed for Titan III exhaust clouds 
(ref- 48) . The most extreme and best documented example out of eight Titan III cases 
observed thus far was for the September 5, 1977, launch (ref. 48). For this launch, 

45 data points for maximum in-cloud HCl concentration C(HCl), obtained over, the 
postlaunch period of 3 to 300 minutes, are well characterized by 3 = 0.54 in the 
empirical expression 


C(HCl) = 55 1 ppmv HCl 


An inclusive envelope is shown in figure 10 which characterizes the previously 
discussed (fig. 8) range of dispersive decay of O and potential rain pH that might 
be expected for Titan III launches. This envelope encompasses all seven standard 
meteorological conditions at Cape Canaveral, but does not include additional uncer- 
tainty bounds to account for the above-mentioned effects. The shaded area beyond 
X = 200 km is meant to denote an additional level of uncertainty, as discussed above. 
The purpose of the envelope is to define a regime of most probable O based on the 
present modeling results. 


Dispersive Decays of In-Cloud HCl Concentration 
and Comparison With Measurements 

Dispersive decays of corresponding vertically averaged HCl concentrations C(HCl) 
were deduced for the eight Titan III cases, and the results are shown as functions 
of X and elapsed time in figures 11 and 12. These averaged peak-centerline HCl 
concentrations were computed from expressions for O by estimating effective cloud 
thickness which applies in figures 4 (a) to (h) for ^ 5 km, as indicated in 

table 1, and by dividing the appropriate expressions for o in table 1 by the corre- 
sponding cloud thickness - 

Since the dispersive decays of C(HCl) in figures 11 and 12 correspond to decays 
of potential rain pH in figures 8 and 9, comparisons of these figure sets facilitate 
preliminary estimates of acid-rain potential based on measurements of in-cloud HCl. 
Such comparisons have proven useful in correlating in-cloud HCl concentrations with 
simultaneous rain composition data. (See ref. 48.) 

A comparison of the calculated results in figure 12 with experimental measure- 
ments of peak in-cloud HCl concentrations for eight Titan III launches is illustrated 
in figure 13 (reproduced from ref. 48) . Since the calculated dispersion results are 
derived from a set of standard meteorologies for the Cape Canaveral area, they have 
no direct relationships to the actual launch cases. Thus, no one-to-one comparisons 
are possible. However, certain first order features of the respective sets are in 
agreement, as discussed below. 

First, to a good approximation all the experimental in-cloud HCl concentration 
data are shown in reference 48 to be adequately characterized by single-term power- 
law decay expressions, such as those shown to apply for the modeled temporal results. 
Next, the total range of measured HCl concentrations is approximately bounded by the 
total range of predicted concentrations after about 0.2 hr (lower limit of model 
validity) up to the indicated termination of each launch data fit (at 0.5 to 5.0 hr). 
Moreover, the respective HCl ranges at various times are large; for example, the 
ratios of highest to the lowest HCl concentration are about 100 for both sets and 
increase after 1.0 hr postlaunch. Despite this apparent intersection of the respec- 
tive experimental and calculated sets of HCl decay histories, four of the measured 
in-cloud peak HCl concentrations decayed significantly slower with time than the 
slowest calculated concentration (3 = 0.83 for FFW) . Notably, the worst-case model 
calculations of HCl concentration (highest values) are partially confirmed experi- 
mentally. The apparent tendency of the dispersion model to overpredict the rate of 
decay of peak HCl concentration is consistent with the derivation of MDM empirical 
dispersion coefficients, which are essentially based on correlations of relatively 
small-scale turbulence measurements. Thus, while there is surprisingly good overall 
agreement with respect to HCl concentrations, the present MDM application appears 
somewhat deficient in that it fails to account for the reduced effect of large-scale 
turbulent diffusion, bounded by an inversion layer, that applies to large stabilized 
SRM exhaust clouds several minutes after launch . 


Reduction of O as a Result of Precipitation Scavenging 


The use of MDM-5(II) to 
HCl(g) washout is checked by 


calculate in-cloud HCl concentrations during progressive 
comparing MDM-5(II) concentration profiles with similar 
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"no rain" results from MDM-4(II) through application of the idealized HCl (g) washout 
model developed in references 6 and 24. The FFW meteorological case is selected, a 
HCl(g) washout coefficient A corresponding to 7.7 mm rainfall per hour is specified 
(A = 4.68 X 10”^ sec“^) , and the elapsed time for onset of rain after cloud stabili- 
zation tj^ is defined (tj_ = 394 sec) . (See ref. 24 for details of the washout 
model.) Two otherwise identical calculations of vertical HCl (g) profiles are needed 
for each downwind distance, one with MDM-5(II) that allows evaluation of O with 
rain and a reference calculation with MDM-4(II) that defines O with no rain. The 
MDM-4(II) gives results equivalent to results presented earlier in figures 5(a) 
to (h) . Equating the ratio of O (with rain) to O (no rain) with the ratio of 
p(HCl) to Pq(HCI) and plotting selected values as a function of distance from 
cloud stabilization X^g gives the results shown in figure 14 (a) . The straight 
line, which agrees very well with the calculated points, was obtained from the washout 
model in reference 24, that is. 


p(HCl)/p^(HCl) 


exp 



^1000 




Note that p (HCl) /p^ (HCl) =1 at = 2.41 km, which corresponds to onset of rain 

at = 394 sec. Results of a similar test of MDM-5(II) using the SFW meteorology 

are shown in figure 14(b) and indicate equally good agreement. The difference in 
downwind HCl (g) washout histories between figures 14(a) and (b) is entirely due to 
the difference in average SRM cloud speed . 

The above comparisons demonstrate that the calculated reduction of O as a 
result of HCl washout, using the numerical outputs of MDM-4(II) and MDM-5(II), is 
quantitatively equivalent to that evaluated with the analytic washout model. However, 
the simplicity of the analytic model, which requires only a single determination of 
a decay in the absence of rain, allows subsequent evaluations of a and various HCl 
washout and deposition characteristics with greatly reduced computational effort. 

Thus, a comprehensive matrix of lengthy MDM-5(II) calculations is no longer required 
in order to evaluate the effects of various combinations of assumed rain onset time 
and washout coefficient corresponding to an assumed rainfall intensity. This tech- 
nique could prove useful for real-time predictions and probability assessments of 
acid-rain hazards during prelaunch countdowns. 


SUMMARY OF RESULTS 

Some important features of the collective results are now summarized. First, 
the dispersive decays of O and C(HCl) differ greatly for the seven standard 
meteorological regimes tested at Cape Canaveral, Florida. A range of more than two 
orders of magnitude in O and C(HCl) is spanned at X ^ 100 km downwind and 
t ^ 2.0 hr. At shorter distances and times, the total span in 0 and C(HCl) is 
somewhat less, but still exceeds an order of magnitude for X > 10 km and t > 0.2 hr. 
These results are consistent with our working hypothesis that various representative 
meteorological conditions at Cape Canaveral lead to widely different exhaust cloud 
dispersion rates. 

Second, significantly large values of 0 and correspondingly large values of 
C(HCl) are shown to be possible out to relatively large downwind distances and times 
from launch site. This is especially true for the two least dispersive Titan III 
cases, which had the largest in-cloud concentrations. Here, for example. 
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0 - 6000 ppmv-m and C(HCl) - 5 ppmv were calculated at X = 50 km and 
t = 5.0 hr. Note that HCl concentrations at ground level are much smaller. Applica- 
tion of a previously developed HCl washout model indicates that the above 0 could 
lead to rainwater pH ^ 1.5/ which is considered environmentally significant for a 
single exposure on certain plants. A pH of 1.5 for rainwater is more than 1000 times 
more acidic than average rain in the Cape Canaveral area (pH = 4.6), which itself is 
10 times more acidic than unpolluted rain. Thus, the 0 values for the two worst- 
case conditions suggest that atmospheric dispersion may not be rapid enough under 
some conditions to avoid significant environmental consequences downwind if rainfall 
occurs . 

Next, the calculated decays of C(HCl) were compared with published analyses of 
available in-cloud peak HCl concentration data obtained from eight Titan III launches. 
Four major features are noted. First, details of the in-cloud data presented else- 
where indicate that straight-line fits of log C(HCl) versus log t are entirely 
adequate representations of the data, which extend beyond t = 2.0 hr for three of 
the eight launches. Thus, the data are in reasonable agreement with the power-law 
characteristics of the present model. Second, the various in-cloud data exhibit a 
spread of two orders of magnitude for t ^ 1.0 hr, depending on launch date and 
meteorology. Thus, the natural variability is large. Third, an inclusive C(HCl) 
envelope bounding the calculated decays of C(HCl) versus t is found to bound 
nearly all the in-cloud HCl concentration data for t ^ 0.2 hr. In fact, the total 
span of in-cloud data at 1.0 hr postlaunch is 0.3 to 20 ppmv HCl (allowing small 
extrapolations) for the eight Titan III launches, which is nearly identical to the 
calculated span for the seven meteorological regimes. Fourth, the power-law decays 
of the in-cloud data tend to be less dispersive than those calculated for the eight 
Titan III cases. More specifically, values of 6 for the in-cloud data vary princi- 
pally from 1.25 to 0.54 for the six best defined cases, whereas values of S for the 
calculated decays of 0 and C(HCl) vary from 1.98 to 0.83 for the seven standard 
meteorological regimes . 

Thus, while some apparent differences exist among the above collective sets of 
experimental in-cloud data and model-calculated values of C(HCl), there is still a 
significant level of overall agreement observed. In effect, this provides a first 
order validation of the model as an assessment tool. However, the present results do 
not allow detailed model-prediction versus experimental-measurement comparisons, 
since the calculated results are based on standard meteorological regimes rather than 
meteorological data for the specific launches. 

Finally, several factors are identified that affect the validity of the calcu- 
lated values of 0 and C(HCl) with respect to the real atmosphere, especially at 
large X and t. First, the model assumes stable stratification conditions in the 
lower troposphere, and hence, lack of vertical convective motion. Obviously this is 
an oversimplification that can significantly invalidate the results under some condi- 
tions. Second, neglect of convective loss of HCl from the SRM exhaust cloud’s upper 
boundary, HCl sorption at ground level, and variable advection for the respective MDM 
layers tends to result in unrealistically large values of 0 and C(HCl) at large 
downwind distances. Neglect of these factors, however, is a convenient means of 
allowing the model to predict somewhat high at large X so as to err on the safe 
side. Finally, the use of variances in the model based on relatively small-scale 
(compared to SRM cloud size) atmospheric turbulence measurements tends to result in 
overestimation of dispersion at large distances and, hence, underestimation of 0 
and C(HCl). The previously cited in-cloud HCl measurement results are considered 
evidence for this effect. 



CONCLUDING REMARKS 


The atmospheric dispersion results for large solid rocket motor exhaust clouds 
calculated with the present multilayer diffusion model indicate wide variability with 
meteorological conditions and sufficiently slow dispersion under some conditions to 
warrant continued environmental concern and further study. In view of these results, 
and since the basic formulation of modified Gaussian dispersion models is statistical- 
analytical in nature and does not treat the essential physics of atmospheric motion 
in the troposphere, it is apparent that a much more comprehensive dynamics model is 
needed. Such a model would provide a more accurate predictive tool for assessment of 
atmospheric transport, dispersion, precipitation scavenging, and dry /wet deposition 
effects. Such a model is especially needed in the vicinity of the Cape Canaveral, 
Florida, area,, where sharp contrasts in surface thermal and moisture characteristics 
tend to control convective activity. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
November 25, 1981 
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APPENDIX 


SUMMARY OF MDM-4(II) INPUT PARAMETERS 

The MDM-4(II) input parameters, used for calculating Titan III dispersion cases 
under each of the seven standard meteorological regimes, are summarized in table Al . 
Note that the present calculations are based on a corrected value for total heat 
release Qj of 2790 cal/g of propellant (see ref. 33) rather than the earlier, 
unrealistic value of 691 cal/g used in reference 30. The cloud centroid height 
shown in parentheses, stemmed from the outdated Qj value and is presently retained 
as an identifier. The unbracketed was calculated directly from the Q(HCl) 

versus (^bk ^TK^/^ data. Note that H, which is defined as the effective source- 
cloud release height for a Model 3 calculation, compares approximately to Other 

Titan III parameters used were a fuel expenditure rate of 4.17 Mg/s (ref. 30); a 
time after ignition t^ in seconds defined as a function of altitude Z (m) by 
tj^ = 0.6346 zO- 4837. a 20.8 weight-percent HCl exhaust composition; and an entrainment 
coefficient of 0.64, used for calculating the width of stabilized Titan III exhaust 
clouds. The Titan III on-pad abort case was also the same as that defined in refer- 
ence 30. The symbols used in table Al are defined as follows: 

H effective source-cloud release height for a Model 3 calculation 

p pressure, millibars (1 millibar = 100 Pa) 

Q(HCl) HCl source strength in each defined SRM exhaust cloud layer, kg HCl 

T temperature, K 

U mean wind speed, m/sec 

UpjQ wind speed at reference height, m/sec 

Z height, m 

Zj^ ground-cloud stabilization height, m 

Zj^ reference height of 2 m 

aj^ lateral diffusion exponent 

vertical diffusion exponent 
0 mean wind direction, deg 

O^ standard deviation of wind azimuth angle, deg 

standard deviation of wind azimuth angle at reference height, deg 
Og standard deviation of wind elevation angle, deg 

standard deviation of wind elevation angle at reference height, deg 
standard deviation of alongwind distribution at source, meters 
UyQ standard deviation of crosswind distribution at source, meters 
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standard deviation of vertical distribution at source ^ meters 


cloud- stabilization time^ sec 

sampling time for standard-deviation calculations, 


$ 


potential temperature 


^ ^1000 


K 


Subscripts : 

BK bottom of Kth layer 

TK top of Kth layer 


sec 
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TABLE Al.- MDM-4(II) MODEL INPUT VALUES FOR STANDARD 
METEOROLOGICAL REGIMES AT KENNEDY SPACE CENTER 
(a) Fall fair weather regime 


Parameter 

Units 

Input values for SRM exhaust cloud layer - 














2 

3 

4 

5 

6 

7 

8 

9 

10 

^BK 

m 

2 

200 

400 

600 

800 

1000 

1200 

1400 

1600 

1800 

^TK 

m 

200 

400 

600 

800 

1000 

1200 

1400 

1600 

1800 

2000 

Zm 

m 





(915) 


1322 




^BK 

m/sec 

3.80 

5.97 

6.39 

6.65 

6.85 

7.00 

6.70 

6.45 

6.18 

5.90 

^TK 

m/sec 

5.97 

6.39 

6.65 

6.85 

7.00 

6.70 

6.45 

6.18 

5.90 

5.60 

CD 

deg 

90 

95,8 

101.6 

107.4 

113.2 

119.0 

121.5 

124.0 

126.0 

129.0 

CD 

deg 

95.8 

101.6 

107.4 

113.2 

119.0 

121.5 

124.0 

126.0 

129.0 

131.0 

^BK 

K 

297.9 

298.2 

298.6 

298.6 

299.2 

299.0 

300.3 

301.6 

302.9 

304.3 


K 

298.2 

298.6 

298.6 

299.2 

299.0 

300.3 

301.6 

302.9 

304.3 

305.6 

^BK 

millibars 

1013 

994 

972 

949 

926 

906 

885 

865 

844 

824 

Ptk 

millibars 

994 

972 

949 

926 

906 

885 

865 

844 

824 

805 

O-p 

^BK 

deg 

10.86 

6.88 

6.42 

6.17 

6.00 

5.86 

,10 

.10 

.10 

.10 


deg 

6.88 

6.42 

6.17 

6.00 

5.86 

.10 

.10 

.10 

.10 

.10 

*^*BK 

deg 

12.00 

7.60 

7.10 

6.82 

6.63 

6.48 

.50 

.50 

.50 

.50 


deg 

7.60 

7.10 

6.82 

6.63 

6.48 

.50 

.50 

.50 

.50 

.50 

Tk 

sec 

364 










''^oK 

sec 

600 










^RK 

m 

18 










Urk 

m/sec 

4.7 










°AR 

deg 

12.0 










^ER 

deg 

6.6 










<^xo 

m 

29.77 

! 89.3 

148.8 

208.4 

267.9 

327.4 

370.7 

311.2 

251.7 

192.1 

®yo 

m 

29.77 

89.3 

148.8 

208.4 

267.9 

327.4 

370.7 

311.2 

251.7 

192.1 

*^zo 

m 

57.16 

57.74 

57.74 

57.74 

57.74 

57.74 

57.74 

57.74 

57.74 

57.74 

“k 


1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

6k 


1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

H 

m 

1273 










Q(HCl) 

kg 

28.6 

1 

101.5 

333.4 

833.2 

1586 

2298 

3116 

2989 

2167 

1423 
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TABLE Al.- Continued 


(b) Spring fair weather regime 


Parameter Units 


Input values for SRM exhaust cloud layer 


m/ sec 
m/sec 


6.95 7.08 

7.08 7.18 


millibars 1013 
millibars 994 


m/sec 


299.4 300.0 300. i 

300.0 300.5 300. ^ 


5.36 4.97 
4.97 4.90 
5.45 5.35 
5.35 5.27 


300.9 301.0 301.0 

301.0 301.0 301.8 


301.8 


29.77 89.30 148.8 208.4 267.9 327.4 387.0 370.1 310.6 251.1 

29.77 89.30 148.8 208.4 267.9 327.4 387.0 370.1 310.6 251.1 

57.16 57.74 57-74 57.74 57.74 57-74 57.74 57.74 57.74 57.74 


84.03 263.5 653.1 1279 


2550 h2m 
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TABLE Al.- Continued 
(c) Low-level sea breeze regime 




Input values for 

SRM exhaust cloud layer 

- 

Parameter 

Units 

1 

2 

3 

4 

5 

6 

7 

^BK 

m 

2 

150 

300 

500 

700 

1000 

1500 

^TK 

m 

150 

300 

500 

700 

1000 

1500 

2000 

Zm 

m 




(567) 

959 



Ubk 

m/sec 

2,51 

7.90 

9.50 

5.60 

4.00 

2.70 

2.9 

Utk 

m/sec 

7.90 

9.50 

5.60 

4.00 

2.70 

2.90 

3.1 

®BK 

deg 

140 

145 

150 

162 

173 

190 

240 

®TK 

deg 

145 

150 

162 

173 

190 

240 

250 

"^BK 

K 

292.2 

293.3 

293.6 

295.5 

297.7 

300.6 

303.2 

^TK 

K 

293.3 

293.6 

295.5 

297.7 

300.6 

303.2 

306.0 

Pbk 

millibars 

1013 

1000 

983 

961 

937 

906 

855 

Ptk 

millibars 

1000 

983 

961 

937 

906 

855 

805 


deg 

9.63 

3.06 

2.54 

.10 

.10 

.10 

.1 

*^Etk 

deg 

3.06 

2.54 

.10 

.10 

.10 

.10 

.1 

Oa 

^BK 

deg 

12.0 

3.81 

3.17 

.50 

.50 

.50 

.5 


deg 

3.81 

3.17 

.50 

.50 

.50 

.50 

.5 


sec 

199 







^oK 

sec 

600 







Zrk 

m 
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^RK 

m/sec 

4.5 







®AR 

deg 

12 







®ER 

deg 

9.9 







^XO 

m 

22.33 

66.98 

119.1 

178.6 

253.0 

168.2 

93.0 

^yo 

m 

22.33 

66.98 

119.1 

178.6 

253.0 

168.2 

93.0 

a 

zo 

m 

42.72 

43.30 

57.74 

57.74 

86.60 

144.3 

144.3 



1 

1 

1 

1 

1 

1 

1 

Sk 


1 

1 

1 

1 

1 

1 

1 

H 

m 

908 







Q(HCl) 

kg 

34.80 

132.5 

545.6 

1587 

3137 

2363 

834.9 
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Parameter 

Units 

^BK 

m 

^TK 

m 

Zm 

m 

%K 

m/sec 

U^K 

m/sec 

®BK 

deg 

CD 

deg 

'^bk 

K 

^TK 

K 

^BK 

millibars ] 

Ptk 

millibars ] 

Ce 

^BK 

deg 


deg 

*^^BK 

deg 


deg 


sec 

^OK 

sec 

Zrk 

m 

Urk 

m/sec 


deg 

^ER 

deg 

^xo 

m 


m 

^zo 

m 

«K 


Bk 


H 

m 1( 

)(HC1) 

kg 




TABLE Al,“ Continued 
(d) Sea breeze regime 

Input values for SRM exhaust cloud layer - 


6.0 

8.9 

150 

150 

289.4 

288.2 


6.70 

4.54 

8.00 

5.42 

247 

600 

2 

6.0 

8.0 

7.6 

14.9 

14.9 

28.3 

1 

1 

)30 

26.2 


8.9 9.6 9.9 10.2 10.4 

9.6 9.9 10.2 10.4 10.6 

150 150 152 153 157 

150 152 153 157 160 

288.2 287.2 286.1 285.2 284.3 

287.2 286.1 285.2 284.3 283.5 

L 004 994 983 972 961 

994 983 972 961 949 

4.54 4.24 4.07 3.95 3.87 


3.95 3.87 

3.87 3.80 


4.24 4.07 3.95 3.87 3.80 
5.42 5.06 4.86 4.72 4.62 
5.06 4.86 4.72 4.62 4.53 


10.6 

10.8 

160 

170 

283.5 

282.8 

949 

937 

3.80 

3.74 

4.53 

4.46 


700 

800 

( 832 ) 

10.8 

10.9 

170 

180 

282.8 

282.2 

937 

926 

3.74 

3.69 

4.46 

4.40 


800 

1300 

993 

10.9 

10.0 

180 

228 

282.2 

287.2 

926 

875 

3.69 

.10 

4.40 

.50 


1300 1800 

1800 2200 

10.0 11. c 

11.9 13 . C 

228 240 

240 250 

287.2 284.2 

284.2 280.2 


44.7 74.4 104,2 134.0 163.7 193.5 223.3 300.6 151.8 93.0 
44.7 74.4 104.2 134.0 163.7 193.5 223.3 300.6 151.8 93.0 
28.9 28.9 28.9 28.9 28.9 28.9 28.9 144.3 144.3 115.5 


46.3 111.2 240.0 1466.5 1815.8 1284 



APPENDIX 


TABLE A1 . - Continued 


(e) Fair weather, pre-cold front regime 


Parameter Units 


Input values for SRP4 exhaust cloud layer 


m/sec 


m/sec 


290.6 295.6 

295.6 297.7 


millibars 

millibars 


m/sec 


(397) 


297.7 

298.1 


298.1 


298.9 299.8 


298.9 299.8 


300.0 


300.0 


300.1 


300.1 


300.7 


18.19 


44.17 


92.0 148.8 208.4 198.6 139.1 93.0 

92.0 148.8 208.4 198.6 139.1 93.0 

52.54 57.74 57.74 79.67 35.80 57.74 


76.98 405.2 1447 2829 13655 9155 h 9Q^ 



APPENDIX 


TABLE Al.- Continued 


(f) Cold front passage regime 


Parameter 

Units 



Input values for SRM exhaust cloud layer ~ 



1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 



^BK 

m 

2 

125 

250 

400 

613 

800 

1000 

1200 

1400 

1600 

1800 

Ztk 

m 

125 

250 

400 

613 

800 

1000 

1200 

1400 

1600 

1800 

2000 

Zm 

m 





(675) 



1230 




«BK 

m/sec 

5.60 

13.0 

15.0 

15.0 

15.0 

13.7 

12.0 

11.3 

10.4 

8.8 

8.0 

^TK 

m/sec 

13.0 

15.0 

15.0 

15.0 

13.7 

12.0 

11.3 

10.4 

8.8 

8.0 

7.0 

®BK 

deg 

41.0 

44.5 

48.0 

49.0 

51.0 

54.0 

59.0 

66.0 

73.5 

80.0 

86.5 

CD 

deg 

44.5 

48.0 

49.0 

51.0 

54.0 

59.0 

66.0 

73.5 

80.0 

86.5 

91.0 

^BK 

K 

295.3 

295.9 

296.9 

296.9 

297.1 

297.7 

298.5 

299.3 

300.0 

300,7 

301.3 

^TK 

K 

295.9 

296.9 

296.0 

297.1 

297.7 

298.5 

299.3 

300.0 

300.7 

301.3 

302.4 

Pbk 

millibars 

1018.6 

1007.5 

990.5 

974.0 

950.0 

929.5 

908.4 

887.5 

866.5 

846.0 

827.0 

Ptk 

millibars 

1007.5 

990.5 

974.0 

950.0 

929.5 

908.4 

887.5 

866.5 

846.0 

827.0 

807.5 

^^BK 

deg 

8.71 

7.64 

7.47 

7.36 

7.26 

7.20 

7.15 

7.11 

7.07 

7.04 

7.02 

^Etk 

deg 

7.64 

7.47 

7.36 

7.26 

7.20 

7.15 

7.11 

7.07 

7.04 

7.02 

6.99 

^BK 

deg 

100.0 

8.77 

8.58 

8.45 

8.34 

8.27 

8.21 

8.16 

8.12 

8.09 

8.06 

^A 

TK 

deg 

8.77 

8.58 

8.45 

8.34 

1 8.27 

1 8.21 

8.16 

! 8.12 

8,09 

8.06 

8.03 


sec 

301 











^OK 

sec 

600 











^RK 

m 

18 











^RK 

m/sec 

8.8 











^AR 

deg 

10.0 











^ER 

deg 

8.8 











^xo 

m 

18.6 

55.8 

96.7 

150.8 

210.3 

267.9 

327.4 

294.4 

234.9 

175.4 

115.8 

^yo 

m 

18.6 

55.8 

96.7 

150.8 

210.3 

267.9 

327.4 

294.4 

234.9 

175.4 

115.8 

^zo 

m 

35.5 

36.08 

43.30 

61.49 

53.98 

57.74 

57.74 

57.74 

57.74 

57.74 

57,74 



1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 



1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

H 

m 

1145 











Q(HCl) 

kg 

25.79 

53.96 

154.5 

483.8 

1184 

2061 ; 

2891 

3301 

2404 

L518 

989.2 


24 



APPENDIX 


TABLE Al.- Concluded 
(g) Post- CO Id front passage regime 


Parameter 

Units 


Input values for SRM exhaust cloud layer 

- 



1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

^BK 

m 

2 

125 

250 

400 

600 

800 

1000 

1200 

1400 

1700 

Ztk 

m 

125 

250 

400 

600 

800 

1000 

1200 

1400 

1700 

2000 

Zm 

m 





(751) 



1341 



Ubk 

m/sec 

4.30 

8.20 

9.00 

9.60 

10.0 

11.0 

11.0 

11.0 

10.4 

8.6 

^TK 

m/sec 

8.20 

9.00 

9.60 

10.0 

11.0 

11.0 

11.0 

10.4 

8.6 

6.0 

CD 

deg 

80.0 

80.5 

82.0 

80.0 

78.0 

75.0 

71.0 

65.0 

57.0 

40.5 

®TK 

deg 

80.5 

82.0 

80.0 

78.0 

75.0 

71.0 

65.0 

57.0 

40.5 

9.0 

^BK 

K 

294.0 

294.7 

295.8 

295.7 

295.2 

295.6 

295.6 

295.6 

295.9 

297.8 

^TK 

K 

294.7 

295.8 

295.7 

295.2 

295.6 

295.6 

295.6 

295.9 

297.8 

303.7 

^BK 

millibars 

1022 

1009 

993.7 

977.0 

954.0 

932.0 

910.6 

890.0 

868.0 

838.0 

Ptk 

millibars 

1009 

993.7 

977.0 

954.0 

932.0 

910.6 

890.0 

868.0 

838.0 

808.0 

^BK 

deg 

8.47 

4.85 

4.42 

4.14 

3.92 

3.77 

3.66 

3.57 

3.50 

.10 

CJp 

Etk 

deg 

4.85 

4.42 

4.14 

3.92 

3.77 

3.66 

3.57 

3.50 

.10 

.10 

^A 

^BK 

deg 

9.00 

5.15 

4.69 

4.41 

4.17 

4.01 

3.89 

3.80 

3.72 

.50 

°Atk 

deg 

5.15 

4.69 

4.41 

4.17 

4.01 

3.89 

3.80 

3.72 

.50 

.50 


sec 

442 










^oK 

sec 

600 










Z'RK 

m 

18 








1 


Urk 

m/sec 

6.0 








1 


^AR 

deg 

9.0 










°ER 

deg 

8.23 










°xo 

m 

18.6 

55.8 

96.7 

148.8 

208.4 

267.9 

327.4 

387.0 

384.0 

294.7 

°yo 

m 

18.6 

55.8 

96.7 

148.8 

208.4 

267.9 

327.4 

387.0 

384.0 

294.7 

^20 

m 

35.51 

36.08 

43.30 

57.74 

57.74 

57.74 

57.74 

57.74 

86.60 

86.60 

“k 


1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

6k 


1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

H 

m 

1420 










Q(HCl) 

kg 

22.46 

34.92 

85.76 

231.1 

569.1 

1125 

1784 

2271 

2999 

2174 

^Q(HCl) 

kg 

1.50 

3.76 

13.28 

56.37 

221.7 

683.7 

1654 

3139 

4917 

5169 

^Zm 

m 







(1132) 


1506 



^With pad abort. 
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TABLE 1.- SRM EXHAUST CLOUD CHARACTERISTICS, AS DEFINED FROM APPLICATION OF MDM-4(II) TO TITAN III 


LAUNCHES, FOR SEVEN STANDARD METEOROLOGICAL CASES AT CAPE CANAVERAL, FLORIDA 


Meteorological regime 

Meteorology 

abbreviation 

Cloud- 
oentroid 
height, 
Z„, m 

Surface 
mixing- layer 
depth, 
m 

HDH Upper 
cloud 
boundary, 
m 

Effective 
cloud thickness 
at X 2 5 km, 
ro 

Mean wind 
speed for 
cloud, 
Uc, nv/s 

Mean wind 
direction 
(North = 0°), 
deg 

Cloud 

stabilization 

time, 

sec 

Calculated 
distance, 
pad to cloud 
stabilization 
(from Oc), 
km 

Empirical 
virtual source 
distance, 

Xo, km 

HCl 

source 
strength, 
mor g 

Fall fair weather 

FFW 

1322 

1000 

2000 

1000 

6.127 

105 

364 

2.23 

10 

14.88 X lo6 

Spring fair weather 

SFW 

1380 

2000 

2000 

2000 

7.0314 

140 

416 

2.93 

4.0 

14.48 

Low-level sea breeze 

LLSB 

959 

300 

2000 

1700 

6.752 

145 

199 

1.34 

3.5 

8.634 

Sea breeze 

SB 

993 

800 

2200 

1400 

9.923 

165 

247 

2.45 

5.0 

10.12 

Pair weather, pre-oold front 

PW, Pre-CP 

875 

200 

1400 

1180 

2.453 

224 

160 

.39 

3.0 

11.89 

Cold front passage 

CPP 

1230 

2000 

2000 

2000 

6.696 

66 

301 

2.02 

4.0 

15.07 

Post-cold front passage 

Post-CFP 

1341 

1400 

2000 

600 

8.712 

69 

442 

3.85 

2.5 

11.30 

Post-cold front passage 
(pad abort) 

Post-CFP 
(pad abort) 

1506 

1400 

2000 

600 

9.264 

69 

342 

3.18 

2.5 

15.86 



TABLE 2.- EMPIRICAL SRM EXHAUST CLOUD PROPERTIES CHARACTERIZING DOWNWIND DISPERSIVE DECAY OF VERTICAL HCl COLUMN DENSITY 


Deduced from application of MDM-4 (II) to Titan III launches for seven standard meteorological cases at Cape Canaveral, 
Florida; the present power- law- decay-modified SFW Shuttle case is shown for both originally used source strength 
m^ = 89 X 10 ^ g HCl (ref. 6 ) and presently reduced source strength m^ = 61 x 10 ^ g HCl consistent with more recent 
^ Shuttle mission design (ref. 55) __ 


Meteorology abbreviation 

— 

Cloud- 

centroid 

height, 

L 

MOM 

upper cloud 
boundary, 
m 

Effective 

cloud 

thickness at 
Xc8 a 5 km, 
m 

Mean wind 
speed for 
cloud, 
Uc, m/8 

HCl 

source 
strength, 
■Hq, 10^ g 

Empirical parameters for vertical HCl column density, <J 

, ppmv-m 

0 » a(Xcs + 

Xq)"® - 

OflT^ 

0 * a(tcs + 

to)"® • 

CQ 

Otr 

ppwv-m 

Xq. 

kiQ 

3 

0, 

ppmv^m 

hr 

3 

1 

Titan III cases 

1 

FFW 

1322 

2000 

1000 

6.127 

14.88 

1.6 X io5 

10 

0.84 

11.9 X io3 

0.453 

0.84 

SFW 

1380 

2000 

2000 

7.0314 

14.48 

2.6 

4.0 

1.98 

.433 

.158 

1.98 

LLSB 

959 

2000 

1700 

6.752 

8.634 

4.3 

3.5 

1.32 

6.35 

.144 

1.32 

SB 

993 

2200 

1400 

9.923 

10.12 

14.0 

5.0 

1.81 

2.17 

.140 

1.81 

FW, Pre-CP 

875 

1400 

1180 

2.453 

11.89 

2.0 

3.0 

.93 

26.3 

.340 

.93 

CPP 

1230 

2000 

2000 

6.696 

15.07 

4.5 

4.0 

1.93 

.968 

,166 

1.93 

Post-CFP 

1341 

2000 

600 

8.712 

11.30 

; .60 

2.5 

1.13 

1.22 

1 .080 

1.13 

Post-CFP (pad abort) 

1506 

2000 

600 

9.284 

j 15.86 

.85 

i 2.5 

1.13 

, 1.65 

.075 

1.13 


Space Shuttle cases 


SFW (original source strength) 



4000 



7.50 

89.0 

14.6 

4.0 

1.64 

6.56 

0.148 

1.64 

iSFW (modified source strength) 

— 

4000 

— 

7.50 

61.0 1 

10.0 

4.0 j 

1.64 

4.49 

.148 

1.64 
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(b) Spring fair weather (SFW) regime (from ref. 51). 


Figure 1.- Continued. 
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(d) Sea breeze (SB) regime (from ref. 30). 
Figure 1.- Continued. 





(e) Fair weather, pre-cold front (FW, Pre-CF) regime (from ref. 51). 


Figure 1.- Continued. 
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(f) Cold front passage (CFP) regime (from ref. 51) . 

Figure 1.- Continued. 
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(g) Post-cold front passage (Post-CFP) regime (from ref. 51). 


Figure 1.- Concluded. 
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(a) Sea breeze (SB) meteorological regime. Surface mixing 
layer depth is 800 m and = 993 m. 


Figure 2.- Dimensions and layer source strengths of the stabilized ground cloud 
of Titan III exhaast products for Qj = 2790 cal/g propellant. 
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(b) Post-cold front passage (Post-CFP) meteorological regime. 
Surface mixing layer depth is 1400 m and = 1341 m. 


Figure 2.- Continued. 
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(c) Cold front passage (CFP) meteorological regime. Surface 
mixing layer depth is 2000 m and = 1230 m. 


Figure 2.- Concluded. 





(a) Space Shuttle and Titan III launches. 

Figure 3-- Summary of some HCl source strength characteristics 
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(a) Fall fair weather (FFW) meteorological regime. 


Figure 4.- Vertical distribution of HCl at various downwind distances from SRM exhaust cloud 
altitude stabilization for Titan III launches under seven standard meteorological 

regimes at Cape Canaveral, Florida. Profiles represent maximum centerline HCl concentra- 
tions calculated from application of MDM-4(II). 


Altitude, 





Figure 




Maximum centerline HCl concentration, ppmv 


(d) Sea breeze (SB) meteorological regime. 


Figure 4.- Continued. 
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(f) Cold front passage (CFP) meteorological regime. 


Figure 4.-' Continued. 
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(g) Post-cold front passage (Post-CFP) meteorological regime with a normal launch 


Figure 4.- Continued. 







Maximum centerline HCl concentration, ppmv 

(h) Post-cold front passage meteorological regime with an abnormal 
(one motor burn on pad) launch (Post-CPP (pad abort) ) . 


Figure 4.- Concluded. 



Vertical HCl column density, (?, ppmv- 



Corrected distance from launch site, + 10, km 

(a) For FFW meteorological regime, == 14.9 x 10^ g HCl 

up to 2.0-km altitude; = 6.127 m/sec . 

Figure 5.- Empirically corrected fits of downwind decays of vertical HCl 
column density calculated for Titan III launches under seven standard 
meteorological conditions at Cape Canaveral ^ Florida. 



Corrected distance from launch site, 

(b) For SFW meteorological regime. m^ = 14.5 x 10^ g 
up to 2.0-km altitude; = 7.031 m/sec. 


Figure 5.- Continued. 





1 10 100 200 
Corrected distance from launch site, X^g + 3-5, km 


(c) For LLSB meteorological regime, m^ = 8.63 x 10® g HCl 
up to 2.0-km altitude; = 6,752 m/sec. 


Figure 5.- Continued, 
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(d) For SB meteorological regime, m^ = 10.1 x 10^ g HC 
up to 2.2-km altitude; = 9.923 m/sec. 


Figure 5.- Continued. 
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(e) For FW, pre-CF meteorological regime, m^ = 11.9 x 10^ g HCl 
up to 1.4-km altitude; = 2.453 m/sec. 


Figure 5.~ Continued. 
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(f) For CFP meteorological regime. m^ = 15.1 x lo^ g H' 
up to 2.0-km altitude; = 6.696 m/sec. 


Figure 5.- Continued. 
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(h) For Post-CFP meteorological regime with an abnormal (one 
motor burn on pad) launch, m^ = 15.9 x 10® g HCl up to 
2.0-km altitude; = 9.284 m/sec. 


Figure 5.- Concluded. 



Distance from cloud stabilization, km 

Figure 6.- Functional relationship used in reference 6 to 
characterize downwind dispersive decay of O for a 
Space Shuttle SRM exhaust cloud under SFW meteorological 
regime. m_ = 89 x 10^ g HCl. 
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Figure 7.“ Virtual-source-distance corrected empirical functions which 
characterize the modified downwind dispersive decay of a for a 
Space Shuttle SRM exhaust cloud under SFW meteorological regime. 
Functions are shown for both the reference 3 (89 x 10^ g HCl) source 
strength/ and the present stabilized-cloud source strength 
(61 X 10^ g HCl) consistent with reference 55. 
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Vertical HCl(g) column density, a, ppmv- 



Distance from launch site, X = + X^, km 


Figure 8.- Dispersive decay of a and potential rain pH as a function of downwind 
distance for eight Titan III cases (solid line) and one Space Shuttle case 
(dashed line) of exhaust cloud dispersion. Results are based on application of 
MDM-4(II) to the seven standard meteorological regimes. (See table 2.) 
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Figure 9.- Dispersive decay of O and potential rain pH as a function of downwind 
drift time for seven Titan III cases (solid line) and one Space Shuttle case 
(dashed line) of exhaust cloud dispersion. Results are based on application of 
MDM~4(II) to the seven standard meteorological regimes. (See table 2.) 
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Figure 10 Envelope for calculated dispersive decay of O and p 
rain pH for all seven Titan III meteorological cases for overl 
transport in the Cape Canaveral, Florida, area. The indicated 
potential rain pH ranges correspond to rainfall rates of up to 
25 mm/hr. 
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Figure 11.- Dispersive decays of volumetric-average HCl concentration calculated 
as functions of downwind distance for Titan III SRM exhaust clouds for seven 
standard meteorological conditions at Cape Canaveral, Florida. 
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Figure 12.- Dispersive decays of volumetric-average HCl concentration calculated 
as functions of downwind drift time for Titan III SRM exhaust clouds for 
seven standard meteorological conditions at Cape Canaveral/ Florida. 
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Figure 13.- Calculated decays (solid line) of in-cloud HCl concentration 
(vertical average/ cloud centerline) compared with empirical fits of 
experimental measurements (dashed line) of peak (^1 sec) in-cloud HCl 
for eight Titan III launches at Cape Canaveral/ Florida. Note that 
the various calculated and experimental HCl decays are not comparable 
on a one-to-one basis / since the respective input and actual launch 
meteorologies differed. (From ref. 48.) 
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Fractional retention of HCl(g) in SRM cloud, p(HCl)/po(HCl) 
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(a) FFW Titan III exhaust cloud. 

Figure 14.- Results of check on MDM-5(II) calculations of 
HCl (g) washout from SRM exhaust clouds. 
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